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A new description for highly nonlinear potential water waves is suggested, where weak 3D effects 
are included as small corrections to exact 2D equations written in conformal variables. Contrary to 
the traditional approach, a small parameter in this theory is not the surface slope, but it is the ratio 
of a typical wave length to a large transversal scale along the second horizontal coordinate. A first- 
order correction for the Hamiltonian functional is calculated, and the corresponding equations of 
motion are derived for steep water waves over an arbitrary inhomogeneous quasi-lD bottom profile. 
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I. INTRODUCTION 

The problem of water waves is one of the classical fields 
of the hydrodynamics, and it has been studied exten- 
sively over many years. Starting from the middle of 90-s, 
in the theory of two-dimensional (2D) potential flows of 
an ideal fluid with a free surface, the so called confor- 
mal variables have been actively employed 0, 0, 0, HI • 
With these variables, highly nonlinear equations of mo- 
tion for planar water waves can be written in an exact 
and compact form containing integral operators diago- 
nal in the Fourier representation. Such integrodifferen- 
tial equations are very suitable for numerical simulation, 
because effective computer programs for the discrete fast 
Fourier transform (FFT) are now available (see, e.g., 
Based on these equations, a significant progress has been 
achieved in the study of nonlinear dynamics of water 
waves, including the mechanism of sudden formation of 
the giant waves 0, H| • Recently, the exact 2D descrip- 
tion has been generalized to the case of a hig hly space- 
and time-inhomogeneous bottom profile [fjllll However, 
the real water waves are never ideally two-dimensional. 
Therefore there is a need of a theory, which could describe 
strongly nonlinear, even breaking waves and, from the 
other hand, it would take into account 3D effects, at least 
as weak corrections. In present work such a highly non- 
linear weakly 3D theory is suggested as an extension of 
the exact 2D theory. It should be emphasized that exist- 
ing approximate nonlinear evolution equations for water 
waves, for example the famous Kadomtsev-Petviashvily 
equation, equations of Boussinesq type ^lj, or the equa- 
tions obtained by Matsuno ^3 1 are valid just for weakly 
nonlinear water waves, but not for overturning or break- 
ing waves. 

It is a well known fact that a very significant difficulty 
in the 3D theory of potential water waves is the general 
impossibility to solve the Laplace equation for the veloc- 
ity potential (p(x, y, q, t), 



<Pxx + <fyy + Vqq = 0, 



(1) 



in the flow region —H(x,q) < y < r}(x,q,t) between the 
(static for simplicity) bottom and a time-dependent free 
surface, with the given boundary conditions 

( f\y= v (x,q,t)=i J (x,q,t), (d(p/dn)\ y =- H ( X}q ) = 0. (2) 

(Here x and q are the horizontal Cartesian coordinates, 
y is the vertical coordinate, while the symbol z will be 
used for the complex combination z = x + iy). There- 
fore a compact expression is absent for the Hamiltonian 
functional of the system, 

n(x,q,t) 



H{r],il)} = ^Jdxdq J (<pl + <p$ + <p%)dy 

-H{x,q) 

+ 9 -J ifdx dq = lC{r,, M + P{ V }, (3) 

(the sum of the kinetic energy of the fluid and the po- 
tential energy in the vertical gravitational field g). The 
Hamiltonian determines the canonical equations of mo- 
tion (see 0, 0, , and references therein) 



m = 



sn 



-V* 



sn 

8rj 



(4) 



in accordance with the variational principle S J Cdt = 0, 
where the Lagrangian is 



C = / iprjt dx dq — Ti. 
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(5) 



In the traditional approach, the problem is partly 
solved by an asymptotic expansion of the kinetic energy 
K, on a small parameter — the slope of the surface (see 
[l3l H^j. and references therein). As the result, a weakly 
nonlinear theory arises, which is not good to describe 
large-amplitude waves (see for a discussion about the 
limits of such theory) . The theory developed in present 
work is based on another small parameter — the ratio 
of a typical length of the waves propagating along the 
x-axis, to a large scale along the transversal horizontal 
direction, denoted by q. Thus, we define e = (l x /l q ) 2 <C 1 
and note: the less this parameter, the less our flow differs 
from a purely 2D flow. The profile y = r)(x,q,t) of the 



2 



free surface, the boundary value of the velocity poten- 
tial tp(x,q,t) = ip(x,T](x,q,t),q,t), and a given bottom 
profile y = —H(x, q) are allowed to depend strongly on 
the coordinate x, while the derivatives over the coordi- 
nate q will be supposed small: \-q q \ ~ e 1 / 2 , \ip q \ ~ e 1 / 2 , 
\H q \~eV*. 

The paper is organized as follows. Sec. II is devoted 
to a general description of the present approach. In Sec. 
Ill, an explicit expression for the first-order correction 
IC^ is obtained, thus we can take into account, in the 
main approximation, weak 3D effects. 



II. GENERAL IDEA OF THE METHOD 

In the same manner as in the exact 2D theory 0, ^| , 
instead of the Cartesian coordinates x and y, we use 
curvilinear conformal coordinates u and v, which make 
the free surface and the bottom effectively flat: 

x+iy = z = z(u+iv, q, t), — oo < u < +00, < V < 1, 

. (6) 

where z(w,q,t) is an analytical on the complex variable 
w = u + iv function without any singularities in the flow 
region < v < 1. Now the bottom corresponds to v = 0, 
while on the free surface v = 1. The boundary value of 
the velocity potential is <p\ v =\ = ip( u , 9; t). In the case of 
a non-horizontal curved bottom, it is convenient to repre- 
sent the conformal mapping z(w, q, t) as a composition of 
two conformal mappings w 1— » £ 1— > z, similarly to works 

z(w, q, t) = Z(((w, q, t),q). (7) 

Here the intermediate function q, t) possesses the 
property Im£(u + 0i,q,t) = 0, thus resulting in the im- 
portant relation 

CO + i, q, t) = £(«, q, t) = (1 + lR)p(u, q, t), (8) 

where p(u, q, t) is a purely real function, and R = i tanh k 
(here k = ~id u ) is the anti-Hcrmitian operator, which 
is diagonal in the Fourier representation: it multiplies 
the Fourier- harmonics pk(q,t) 
Rk = itanhfc, so that 



/ p(u, q, t)e tku du by 



dk 

Rp(u,q,t) = [itanhk]p k (q,t)e lku — 

Zn 

p(u, q, t) du 



= P.V. / , g x 

J 2sinh[(7r/2)(t2- U )] 1 ' 

(P.V. means the principal value integral.) A known an- 
alytical function Z(£ 7 q) determines parametrically the 
static bottom profile: 

XW(r,q)+iYW(r,q) = Z(r,q), (10) 

where r is a real parameter running from —00 to +00. 
The profile of the free surface is now given (in the para- 
metric form as well) by the formula 

X [s] (u,q,t) + iY [s] (u,q,t) = Z [s] (u,q,t) = Z(£(u,q,t),q). 

(11) 



For equations to be more short, below we do not indicate 
the arguments (u,q,t) of the functions ip, £ £ (the over- 
line denotes the complex conjugate). Also, we introduce 
the notation Z'(£) = d^Z(^,q). The Lagrangian of the 
system in the variables tp, £, and £ can be re-written as 
follows (compare to 9]): 



£ = / Z'(0Z'(0 



£,t£,u — £t 



21 



ip du dq 



ic{ip,z(o,z(m 



21 



A 



2i 



R 



2 

du dq, 



du dq 
(12) 



where the indefinite real Lagrangian multiplier A(u,q,t) 
has been introduced in order to take into account the 
relation (JSJ . Equations of motion follow from the varia- 
tional principle 5 A — 0, with the action A = J Ldt. So, 
the variation by 5tf> gives us the first equation of motion 
— the kinematic condition on the free surface: 



(13) 



Let us divide this equation by Z'(£)| 2 |£ti| 2 an d use the 
analytical properties of the function £t/£ M . As the result, 
we obtain the time-derivative-resolved equation 



Zt = UT + i) 



(SK/Sip) 



toi 2 i& 



(14) 



where the linear operator T = R^ 1 = —i coth k has been 
introduced. Further, the variation of the action A by St; 
gives us the second equation of motion: 



2i 



-i-Im(z(0)\Z'(0\% 



(1 + £JR)A 



2i V °V Wl ™ 2i 
After multiplying Ea. (|15() by — 2i£ u we have 

{[i>t + gi^z(C]\iu\ 2 - ^uitZu) \z'(0\ 2 
'SIC' 



(15) 



(l + iK)A 



2/ 1 ^ ) Z'iOHu, 



(16) 



where A is another real function. Taking the imaginary 
part of Ea. (|16f) and using Ea. (|13|) . we find A: 



A = T 



SIC 

Sip 



2TRe 



SJC 

sz 



z'm, 



(17) 



After that, the real part of Ea. (|16|) gives us the Bernoulli 
equation in a general form: 



tpt + gImZ(£) = tf) u T 
2 Re 



(SlC/Sip) 



\z'(0\ 2 \U 



t[msic/s^)} 
\z'(0\ 2 \Zu\ 2 



(ft - i) {{S1C/SZ)Z I {^ U \ 

WW\U 2 



.(18) 
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Equations l|14|) and (|18fl completely determine the evo- 
lution of the system, provided the kinetic energy func- 
tional lC{ip, Z, Z} is explicitly given. It should be empha- 
sized that in our description a general expression for /C re- 
mains unknown. However, under the conditions \z q \ <C 1, 
\tp q \ <C 1, the potential </?(«, v, q, t) is efficiently expanded 
into a series on the powers of the small parameter e: 



<p = <pW +( pW +tp W + 



(n) 



(19) 



where i^™"*" 1 ) can be calculated from pW, and the 
zeroth-order term ip^ = Re(j>(w, q,t) is the real part of 
an easily represented (in integral form) analytical func- 
tion with the boundary conditions Re</>| u= i = ip(u,q,t), 
Im0|i, = o = 0. Correspondingly, the kinetic energy func- 
tional will be written in the form 

K = /C (0) + /C (1) + /C (2) + . . . , /C (n) -e™, (20) 

where K.(°'{ip} is the kinetic energy of a purely 2D flow, 

IC (Q) W = \ f[(^) 2 + (^) 2 ]dudvdq 



ipRipu du dq, 



(21) 



and other terms are corrections due to gradients along q. 
Now we are going to calculate the first-order correction 



III. THE FIRST-ORDER CORRECTIONS 

As the result of the conformal change of two variables, 
the kinetic energy functional is determined by the expres- 
sion 

^ = \jWl + ^l + J{Cl- V^) 2 ] du dv dq, (22) 

where the conditions x u = y v , x v — —y u have been taken 
into account, and the following notations are used: 

J=\z u \ 2 , (Q • V<p) = aip u + btp v + ip q , 



Consequently, the Laplace equation in the new coordi- 
nates takes the form 

V • (QJ(Q ■ V^)) - 0, (23) 
with the boundary conditions 

<p\ v =i = ip(u, q, t), [cp v + bJ (ifi q + aip u + b(p v )]\ v = = 0. 

(24) 

Under the condition e <C 1 it is possible to write the 
solution as the series (|19|) . with the zeroth-order term 
satisfying the 2D Laplace equation 

fu2 + fvv = 0, <p\v=i = i/>(u,q,t), <Pv\v=0 = 0- 



Thus, it can be represented as = Re <p(w, q, t), where 

Mq,ty kw dk 



4>{w,q,t) 



cosh k 27r ' 



(25) 
(26) 



Ml,t) = J ^(u,q,t)e~ tku du. 

On the free surface 

(j)(u + i, q, t) = 9{u, q, t) = {l + iR)ip(u, q, t). (27) 
For all the other terms in Ea. (|19H we have the relations 

ft +1) + + V • (QJ(Q • V n) )) = (28) 

and the boundary conditions ( ^( T1 + 1 ) |^ =1 = 0, 

+ fc/foW + mp«> + M, n) )]|„=„ = 0. 

Noting that f(tpu tpu + fv (fv)dudv = (it is eas- 
ily seen after integration by parts), we have in the first 
approximation 

JC {1) = \J Jiyf+a^+byWfdudvdq 

du dv dq. (29) 



Re ^ 


, 4>uZq\ 







2 " 



Since z(w) and <fi(w) are represented as z{u + iv) = 

e Hi-v)Zl°\( u ) and <j>(u + iv) = e^ 1 ""^^), we can use 
for u-integration the following formulas: 



du I < 

o 

2k _ ! 



k ( 1 -"U(u)][e*( 1 - t ')B(tt)] dv 
dk 



2k 



A k B k 



2tt 



= - l -J B(u)B?A(u)du 
+ l - ( Wfajd^AW^du, 



(30) 



with A^(u) = e k A(u), BW(u) = e k B(u). As the re- 
sult, we obtain from Ea. l29f) the expression of the form 

KM = /Cg - /Cohere K$ = ^{*, Z, Z}, fc[J> = 
T{^ b \W\,Z^\W\}, with Z = ZW, ZW = e k Z, 
vj/I 6 ] = e % = [coshfc] _ V- The functional T is defined 
below: 



T = - / {Z u y q -Z q y u )d- L (Z u V q -Z q Vu)dudq 



- I { [(Z u Vq - Z q * u f/Z u \ z 



-Z [(Z u V q - Z q y u )yZ u ] du dq. (31) 
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From here one can express the variational derivatives 
(SfC^/Sip) and (5K.W/5Z) by the formulas 



8KP 



Si/j 



5Z 



[cosh fcl 



5T 



77 + (1 - 




( 5JC {1) 









,-fe 



SIC 



(1)' 



(32) 



(33) 



The derivatives {5T/5Z), [8K$ /8&% and 

(SK.^ /SZ^) are calculated in a standard manner, for 



instance, 
6* ~ 8 9 



(Z u * 9 - Z 9 *. u ) + - Z q V u /Z u )Z] 



*™ ft. 



(Z u *, - Z,* u ) + d u [(* g - Z q V u /Z u )Z] 

a- 1 (z u * g -z g * u ) + (* g - z g *„/z„)z 



2 

16 



a u [(# g - z 9 *„/z u ) 2 z] - (* g - z g *„/z u ) 2 z u 



Now one can substitute (SK/Stp) » -.R^u + (SK^/Sip) 
and (8K./8Z) w (SJC^ /SZ) into the equations of motion 
and JTSJ, keeping in mind that Z = Z(£,q), Z u = 
Z'mu, Z q = Z\Qi q + d q Z, = Z([coshk]-y,q), 
and so on. Thus, the weakly 3D equations of motion are 
completely derived, and our main goal is achieved. 



The answers are more compact in the limit |fe| ^> 1, 
corresponding to the "deep water" , when R — > H, T — » 
—H, with H being the Hilbert operator: H — i sign k. In 



this case /C 



(i) 



[b] 



K. 



deep 



0, and therefore 
^H^ u dudq + F{<S!^. i Z,~Z}. (34) 

After appropriate rescaling of the variable u, we can write 

Z = «+(*- H)Y(u,q,t), Z u = l + (i-H)Y u . 

The equations of motion for quasi-plane waves on the 
deep water look as follows: 



Z t = Z U {H - i) [Hi/> u ~ (6F/6*P)]/\Z u f 



(35) 



^ t + gY = ip u H 



[Hip u - (6F/6i/>)]/\Z u 



H 



2Re((H + i)[Z u (ST/5Z)])/\Z u \ 2 1 (36) 



where (SF/Sip) = 2 Re (1 - iH){5T/5^) 
IV. SUMMARY 



Thus, now we have nonlinear evolution equations for 
weakly 3D steep water waves, as for deep water case, 
as for arbitrary quasi-lD bottom profile. The obtained 
equations are intended to describe, for example, the sud- 
den formation of giant waves in open sea, as well as over- 
turning waves on a beach. The following step should be 
development of an efficient numerical method for simu- 
lation these equations. 
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